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Abstract 

We study the role of QCD plasma instabilities in non-equilibrium quark-gluon plasmas. First, 
we argue that such instabilities must drastically modify the "bottom-up" thermalization scenario 
for heavy-ion collisions. Second, we discuss conditions for the existence of instabilities in a more 
general context than previously treated in the QCD literature. We also give a thorough qualitative 
review of the origin of instabilities. We discuss some mechanisms whereby the growth of plasma 
instabilities saturates. Finally, we solve explicitly for instabilities and their growth rates for two 
extreme cases of anisotropic non-equilibrium plasmas that can be treated relatively simply and 
analytically: f{p) = FiPl.) ^{Pz) and /(p) = FiPz) S^'^\p±), where f{p) is the distribution of 
particles in momentum space. 



I. INTRODUCTION 



Mrowczynski [1-6] has long advocated that plasma instabihties may play an important 
role in the thermalization of quark-gluon plasmas from anisotropic, non-equilibrium initial 
conditions, such as those pertaining to heavy-ion collisions. (For other early discussions, 
see Refs. [7-11].) Since heavy-ion collisions are rather complicated, with a large variety 
of different physical processes playing roles at different stages of the collision, it is useful 
to examine this proposal in a theoretically clean limit. A natural framework for such an 
investigation has been provided by Baier, Mueller, Schiff, and Son [12]. In a beautiful 
paper, they consider what happens shortly after the collision in the saturation scenario [13- 
16], at energies high enough that the effective strong coupling a at the relevant momentum 
scale Qs (known as the saturation scale) is weak, a{Qs) <^ 1. The result of their analysis, 
which found an interesting variety of parametrically different time scales characteristic of 
different stages of thermalization (such as a~^l'^Q~^ ^ or^l'^Q~^ ^ and or^'^I^Q~^\ is known as 
bottom-up thermalization. We will show that even the very first moments of this scenario 
are drastically modified by the effects of plasma instabilities, which were not considered by 
Baier et al. 

We will use coUisionless kinetic theory coupled to soft classical gauge fields (also known 
as the Hard Thermal Loop approximation) to study the plasma instabilities during the first 
phase of bottom-up thermalization. First, in section I A, we will briefly review the first 
phase of the bottom-up thermalization scenario and summarize our argument that it must 
be drastically modified. We also argue that the coUisionless approximation is appropriate to 
check for instabilities. For those readers more interested in the plasma instabilities of kinetic 
theory than in the details of the bottom-up thermalization scenario, the upshot is that we 
will be motivated to study the existence and properties of plasma instabilities for an initial, 
homogeneous, non-equilibrium, phase-space distribution of hard gluons of approximately 
the following form: 

fip,x)^F{p^)S{p,), (1.1) 

where p± = Y^p|~-Fp|. This is an interesting limiting case of certain classes of distributions 
that have been studied in the literature numerically by Randrup and Mrowczynski [6] and 
Romatschke and Strickland [17], but (1.1) has the advantage of being more analytically 
tractable.^ The particular plasma instabilities that will be of interest are known in various 
contexts as (i) Weibel [19], filamentation, or pinching instabilities, and (ii) two-stream or 
Buneman [28] instabilities, with the names somewhat interchangeable in each category. Be- 
cause many readers will not be intimately familiar with these instabilities in kinetic theory, 
we will review them both formally and qualitatively in section II for parity-invariant situa- 
tions. Typically, for parity-invariant particle momentum distributions, plasma instabilities 
exist (in the coUisionless approximation) whenever those distributions are anisotropic. We 
will give conditions which make this statement precise. We also discuss certain aspects of the 
physics that eventually cuts off the growth of such instabilities and point out the qualitative 
differences between the cases of Abelian and non-Abelian gauge theories. To make the rather 



Romatschke and Strickland [17] do obtain analytic results in a more general class of situations, but their 
results are unwieldy enough that they refrain from giving them explicitly in their paper. Also, another 
example of analysis of instabilities in ultra-relativistic plasmas, generalized to inclusion of a background 
magnetic field but restricted to wave vectors in certain symmetry directions, may be found in Ref. [18]. 
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FIG. 1: The standard picture of lines of constant proper time r = — in a high-energy 
coUision, in which the 45° lines are the colliding nuclei. The two arrows show how free-streaming 
partons leaving from the same initial point would later segregate themselves in z according to their 
value of Vz- 

general discussion in section II more concrete, we will then analyze in detail in section III 
the instabilities associated with the particle distribution (1.1). This distribution is simple 
enough that much of the analysis can be done analytically. We also give a specific example 
in section IV of how non-Abelian interactions can cut off the growth of plasma instabilities. 
Various matters are left for appendices, including an analytic analysis of instabilities for the 
distribution 

= (1.2) 

which is complementary to that of (1.1). 

We will show that plasma instabilities play a crucial role in the post-saturation evolution 
of heavy-ion collisions, in the framework originally established for bottom-up thermalization, 
but we are not yet able to give detailed parametric estimates of the time scale for equilibra- 
tion (analogous to the a'^^I^Q'^ time scale of Baier et al.). We do expect thermalization 
to be more rapid, but this is a topic for future research. 

A. The first phase of bottom-up thermalization 

Suppose the initial nuclei are moving in the -^z directions. The starting point for the 
bottom-up thermalization scenario is to consider an initial distribution of particles that is 
(i) boost-invariant in the z direction, (ii) homogeneous in the x and y directions, and (iii) 
dominated by gluons with momenta p of order some scale Qs and with non-perturbative 
occupation numbers per mode of order 1/a. For sufficiently high energy collisions, the ffist 
assumption should be adequate for describing the central rapidity region and the second 
assumption for describing local thermalization for large nuclei. These two assumptions 
imply that physical quantities (measured in a local rest frame) depend only on the proper 
time r = \Jt'^ — z'^, as depicted in Fig. 1. The third assumption is the saturation assumption 
and is taken to hold at r ~ Qs^- The bottom-up thermalization scenario concerns itself 
with what happens at later times r ^ Qs^, given this initial condition. Following Baier et 
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ai, we will refer to the initial ghions as "hard" gluons.^ 

If there were no interactions, then, as time progressed, the hard gluons would simply 
follow straight-line trajectories and segregate themselves in z according to their values of 
''^z = Pz/p, as depicted in Fig. 1. One can show that, in a local rest frame (where the average 
p vanishes), the resulting distribution of hard gluons for r ^ would then be dominated 
by gluons with Px ^ Py ^ Qs but ~ 1/r <^ Qg. The occupation numbers of these modes 
would still be of order 1/a. The spatial number density A^h of hard gluons decreases linearly 
with the linear expansion as 

Baier et al. argue that small-angle scattering between the hard gluons instead widens the 
spread in to 

~ '^-'^ 

during the first stage, 1 -C QsT -C a'^^"^, of their scenario. They then argue that hard gluon 
number is conserved over the short duration of this first stage, and so the spatial number 
density of hard gluons still decreases linearly with r as (1.3). The occupation number per 
mode (related by — /h ~ PxPyPzfh) then decreases as 

QIPz «(ar)V3- ^'■^) 

The fact that /h <^ l/a for QgT 3> 1 means that the direct interactions between hard 
particles can be treated perturbatively at these times. 

Unscreened Coulomb scattering is infrared divergent. The small-angle scattering of hard 
gluons discussed above is therefore dominated by momentum exchanges of order the inverse 
screening length. Baier et al. estimated this soft scale, which we will call mgoft, to be^ 



^soft 



P Qs r ^ ' 



They loosely referred to mgoft as the Debye mass mD and took estimates from the literature 
[20-22]. In fact, the Debye mass is not precisely the correct quantity to use since one really 
wants to use the full self-energy 11'''' (a;, q) for the soft exchanged gluon of the small-angle 
scattering process. Also, the Debye mass is not precisely defined in anisotropic situations: it 
depends on direction. However, these details do not affect parametric estimates, and so (1.6) 
is a perfectly adequate estimate of the relevant soft physics scale. There is only one problem. 
A parametric estimate does not necessarily reveal the sign of the result. As we shall see, 
a more appropriate definition of ml^^^ reveals that its sign is negative for certain modes. 



^ This usage is to distinguish them from lower energy partons generated later in the evolution of the collision. 
Prom the point of view of the colored glass condensate picture of the nucleus, or of "mini-jet" analyses 
of heavy ion collisions, the usage may seem peculiar, since Qs is treated as an infrared cutoff in those 

contexts. 

^ Our terminology is slightly different from that of Baier et al. [12]. They refer to (1.4) as the "soft" physics 
scale during the first stage of their scenario. In our paper, the term "soft" will instead generally refer to 
(1.6). 
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corresponding to an instability in the soft gauge fields in the presence of the anisotropic 

distribution of hard gluons. These unstable soft modes then grow exponentially with a 

characteristic time scale of 

1 [T' 

|"^soft| V 

Because the growth time tgrowth is parametrically small compared to the expansion time 
T for QgT ^ 1, the first stage of bottom-up thermalization is drastically modified by the 
generation of large, non-perturbative soft gauge fields. Since this is true for any QgT ^ 1, 
the first stage of the original bottom-up thermalization scenario cannot, even briefly, be 
qualitatively accurate. 

In order to demonstrate the instability, note that the first stage of the bottom-up scenario 
is described by hard gluons with Pz Px ^ Py ^ Qs as in (1.4). Using the azimuthal 
symmetry of the problem, such a distribution of hard gluons can be approximately described 
as having the form (1.1). In Sec. Ill, we will analyze the soft instabilities created by hard 
distributions of this form. 

The extreme planarity (oblateness) of the anisotropy of our model distribution (1.1) is 
somewhat an oversimplification of the actual situation in the bottom-up scenario. Baier et 
al. point out that not only hard gluons contribute to rri^^f^, but so also do softer gluons which 
have been generated by Bremsstrahlung and then pushed to momenta kg of order pz (1.4) 
by scattering from hard particles. These softer gluons will be anisotropic but not have the 
extreme anisotropy of the hard gluons. In the bottom-up thermalization scenario, Baier et 
al. show that such soft gluons give a contribution to ml^f^ comparable to that of the hard 
gluons. We will see that whenever there is significant deviation from anisotropy, then there 
are unstable modes with momenta of order mgoft, and so plasma instabilities will indeed 
modify the bottom-up scenario. However, as a simple analytically-tractable example of this 
phenomena, it will be instructive to ignore the softer contribution and study the effects of 
just the extreme oblate distribution (1.1) of the hard gluons. 

In our analysis of plasma instabilities and their growth, we will use a coUisionless approx- 
imation for the hard particles. We need to justify this approximation since a significant rate 
of collisions could eliminate the instability. Baier et al. estimate the rate of hard particle 
collisions for a single hard gluon in the first stage of their scenario as 

— ~(7A^h(l + /h) 2 ' (1-8) 

where a ~ a^ml^^^^ is the cross-section. Using the previous estimates for N-^, msott, and Pz, 
the mean free time between such collisions is then 

For QsT ^ 1, this is parametrically long compared to the time scale tgrowth of (1-7) for 
the growth of instabilities. We may therefore ignore collisions of hard gluons for the pur- 
pose of understanding whether the bottom-up thermalization scenario is affected by plasma 
instabilities. 

In our analysis of instabilities, we will assume that the plasma is homogeneous. In fact, 
from Fig. 1, one can see that densities will vary in a constant t slice over a distance scale of 
order r. The homogeneity assumption is valid for QgT ^ 1 because this scale will be large 
compared to the typical wavelength 1/| mgoft | of unstable modes. 
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II. PLASMA INSTABILITIES 



The basic starting point is the expression for the soft gauge-field self-energy 11^'' (a;, q) in 
the presence of a given homogeneous but non-equilibrium distribution /(p) of hard particles. 
In this section, we will treat general distributions f{p) and will not yet specialize to specific 
forms such as the planar distribution of (1.1). 



A. Review of the retarded self-energy 

First, let us briefly review the known result for the self-energy from coUisionless kinetic 
theory. In order to retain contact with the literature on non-relativistic plasmas, we will 
avoid specializing to the ultra-relativistic limit for the moment. For simplicity, we start by 
discussing Abelian gauge theories such as electromagnetism. One starts with the Vlasov 
equations, which are the coUisionless Boltzmann equation for a hard particle distribution 
f{p,x, t) and Maxwell's equations for the soft gauge fields: 

dtf + v- VJ + e{E + vxB)- V^/ = 0, (2.1) 
Q^F'^-^ji^^e f v'^f, (2.2) 

where e{E + v x B) is the Lorentz force on hard particles of charge e, v = v{p) is the hard 
particle velocity for a given momentum p, and v'^ = {l,v). There is also an implicit sum over 
particle species and their charges on the right-hand side of (2.2). In the absence of external 
fields, one treats E{x,t) and B{x,t) as small and hncarizcs f{p,x,t) = fo{p) + fi{p,x,t) 
in small fluctuations /i about some homogeneous distribution fo{p)- We will assume that 
e/pW^/o vanishes (when implicitly summed over species) so that /o carries no net charge or 
net current. Then 

dtfi + V ■ V,/i + eiE + vxB)- Vpfo ^ 0, (2.3) 

a.F'^^ = e [ v^h, (2.4) 
Jp 

which are known as linearized Vlasov equations. Fourier transforming from (i, x) to (a;, q), 
one may solve for /i and obtain 

,e„^..,,e=/l^(£±^ii£i^, (2.5) 
Jp -u + v -q-ie 

where the e is a positive infinitesimal inserted as a prescription to obtain the retarded 
solution. We use the shorthand notation 

'''' ■ ■ (2.6) 



p 



(27r) 



We may rewrite (2.5) as 

iQ.F^''' ~ -n'^M,, (2.7) 
extracting the result for the retarded self-energy 11'^'^, 

dfip) 



Jp 



dp'' 



-v^'g'"" H 

—CO + v-q — i€ 



(2.8) 
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where we use ( —+++) metric convention, and Roman indices such as k are spatial indices 
that run from 1 to 3. The self-energy U^'^ is symmetric in the indices fiu^ It is also useful to 
note that the above expression satisfies the Ward identity, QnH^^ = (which is just current 
conservation, Q^j^ = 0), where Q = {io,q). This implies 



^o.^9_LJ_^ noo = ^-V-- (2-9) 



Id *^citn*^n-J-i 11 T r i j-1 ^-^n -J- 1 r» 1 r\'^r^r\t-\-i- -J-l-irv-J- j-i^ 



The result for non-Abelian gauge theory^ is essentially identical except that e is replaced 
by g"^ tr(T"Tj') = S'^^Cgds/ for each species s of hard particle, where g is the gauge coupling, 
a and h are adjoint color indices for the soft gluon, are the color generators for the 
hard particle's color representation, Cs is the quadratic Casimir given by T"T" = C^, and 
dsi dx are the dimensions of the particle's representation and the adjoint representation 
respectively. For gluons in QCD, Cgds/dx — 3, while for quarks Csds/d^ — 1/2. However, 
since there is no essential difference between the Abelian and non-Abelian results of the 
linearized analysis, we will generally write the self-energy as (2.8) for the sake of brevity, 
with the understanding that 

e^-E-^^ (2.10) 

s 

with f ^ fs, where represents the number of degrees of freedom of a given hard particle 
species, excluding color. (So, for instance, if /g for gluons represents the density of gluons 
per spin state and color, then Ug would be 2.) 

Isotropic distributions f{p) = f{p) have the following special properties: (i) 11(0, q) is 
diagonal with the only non-zero entry ll^^{0,q) = —m^, where m-D is the Debye screening 
mass (inverse Debye screening length); and (ii) n(a;,0) is diagonal with the only non-zero 
entries n'''(u;,0) = rri^^S^^, where rripi = niD/^/S is the plasmon mass (plasma frequency). 
In contrast, for generic anisotropic distributions, (i) 11(0, q) is not diagonal and depends on 
the direction q, and (ii) n(a;,0) is not diagonal and the eigenvalues of its spatial part are 
not all the same. 



B. Criteria for Instability 

In momentum space Q — {u!,q), the linearized effective equation (2.7) for the soft gauge 
fields is 

[(_^2 ^ ^2^g^. _ Q^.QU ^ g)]^^ _ Q_ (2.11) 

For a given q, there will be instabilities if there are any solutions for u with Imcu > 0. 

For the case of parity-invariant background particle distributions, f{p) = f{—p), there 
are some simple sufficient conditions for the existence of instabilities. In the rest of this paper, 
we will assume that distributions are parity invariant.^ There is then a simple sufficient (but 
not necessary) condition for the existence of an instability: 



The symmetry is not obvious for the first term in (2.8) unless one integrates by parts and uses v = 'VpEp 
so that dpiV^ = dpi dpi Ep, where Bp is the energy of a hard particle with momentum p. 

^ See Ref. [5] for the non-Abelian version of (2.8) in the non-equilibrium case. However, the basic modifi- 
cations have been known in various approximations for ages {e.g. Ref. [23]). Ref. [24] also has expressions 
for (2.8) specialized to the case of axial symmetry. 

^ We also implicitly assume that hard particle dispersion relations are parity invariant. 
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FIG. 2: Example of an effective potential V^ff (^) iii scalar field theory for which $ = would be 
unstable. 



Condition 1. There is an instability associated with a given wavenumber q for 
each negative eigenvalue of the 3x3 matrix q^5'^^ + n*-' (0, q). 

This is a simple generalization of the Penrose criteria used by Mrowczyhski in Ref. [2] for 
a special case. There is a simple continuity argument for this assertion, which we give in 
Appendix A 1.^ 

The zero-frequency self-energy 11(0, q) depends only on the direction q and not the mag- 
nitude q ol q because (2.8) can be rewritten for a; = as^ 

n-(o,,) = e=/^Uv + ^ 

Jp dp'' I v q- 

We will write n(0, q) instead of n(0, q) to emphasize this. A simple corollary of Condition 1 
is then that there will always be an instability in some mode {i.e. for small enough q) if 
n*-' (0, q) has a negative eigenvalue. 

Readers may find it useful to consider as an analogy the static effective potential Ves{^) 
for the softest modes of a scalar field theory. If T4ff($) has negative curvature at $ = 0, as 
depicted in Fig. 2, then $ = is unstable. For scalars with no intrinsic mass (analogous to 
gauge bosons), this curvature is simply 11(0). The condition 11(0) < implies that there is 
an instability, even though the precise calculation of the growth rate would require studying 



(2.12) 



Ref. [2] considered the special case where q points along a symmetry axis of /(p) and the eigen- modes 
of n'-' are determined by symmetry. Textbook discussions of the Penrose criteria for similar special cases 
in the non-relativistic context may be found in Sec. 9.10.2 of Ref. [25] and Sec. 3.3.4 of Ref. [26]. We 

are unaware of a general analysis of the case where q does not point in a symmetry direction of the 
problem, and so we give our own argument in the appendix. It also sets up the argument for Condition 
2 below. Similarly, one often sees the Penrose criteria playing the role of necessary (as well as sufficient) 
conditions for instability in certain situations, such as where q is in the z direction and f{p) has the 
form Fz{pz)F±{p±) with F{pz) a monotonic decreasing function of \pz\ [2, 25, 26]. We do not know how 
to prove necessity in the case where q points in non-symmetry directions, for reasons described in the 
appendix. 

Physically, the independence of 11(0, q) on the magnitude q of the soft field wave number depends on 
q being small compared to hard particle momenta — an approximation implicit in our use of the Vlasov 
equations. 
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the dynamics, depending on Linearizing ||V$p + Kfj($) about $ = 0, one would 

find that the condition for a particular mode to be unstable is + 11(0, q) < 0, analogous 
to Condition 1. The picture of Fig. 2 also suggests that non-linearities in the interactions of 
the soft fields may eventually cut off the growth of the instability, a topic we will return to 
later. 

We will refer to instabilities discovered by Condition 1 as "magnetic" instabilities because 
the condition itself involves the self-energy for A{u!=0) and so involves magnetic but not 
electric fields. We will give a physical picture of the origin of the instability in Sec. II C. 
Of course, the growth of magnetic fields implies the generation of electric fields as well, by 
Maxwell's equations, so the actual growing mode will not be purely magnetic. 

Note that the Ward identity q^Tl^^ = implies q^W^{0,q) = 0. Therefore, the longitu- 
dinal polarization of A(a;=0) is always associated with a zero eigenvalue of n*-^ (0, q) and so 
will not play a role in Condition 1. We could have phrased Condition 1 equivalently using 
its transverse projection: 

Condition 1 '. There is an instability associated with a given wavenumber q for 
each negative eigenvalue of the 3x3 matrix q^5'^^ — q^q^ + 11*^(0, q). 

In Appendix B, we expand upon an argument in Ref. [27] to use Condition 1 to show 
that some unstable mode always exists unless the matrix W^i^O.q) is identically zero. We 
further show that this leads to a simple necessary as well as sufficient condition for magnetic 
instability in the ultra-relativistic limit: 

Condition 1-b. Magnetic instabilities exist for a given (parity symmetric) distri- 
bution f(p) if 



is anisotropic in p. Moreover, in the ultra-relativistic limit, this is a necessary 
(as well as sufficient) condition for the existence of instabilities. 

We note in passing that the angular average of A4 (p) over p gives 



which defines a basic mass scale we will encounter repeatedly. In the ultra-relativistic case, 
it also turns out to be the effective mass of transverse electromagnetic waves in the large q 
limit, with o;^ ~ -|- m^. This is shown in Ref. [27] and inspires the notation m^, though 
this correspondence only holds in the ultra-relativistic case. Note that is independent 
of the direction q. 

In Condition 1-b we have re-emphasized the underlying assumption of our discussion that 
f{p) is parity symmetric. The condition is obviously incorrect without this assumption: 
Consider, for example, a standard equilibrium distribution boosted to an inertial frame 
where the average momentum does not vanish. The corresponding f{p) and Ai{p) are then 
anisotropic and also not parity invariant, but equilibrium is stable, regardless of the inertial 
frame. 

There is another type of instability that can occur, which can be associated with the 
electric potential at ui = [which for non-zero u is related to the longitudinal polarization 
of A by the Ward identity (2.9)]. In Appendix A 2, we show that another (independent) 
sufficient condition for the existence of an instability is 




(2.13) 




(2.14) 
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Condition 2. There is an instability associated with a given wavenumber q if 
- n°°(0, q) + n°^(0, q) [q' + n(0, q)] n^°(0, q) < 0, (2.15) 

where [g^ + 11] ~^ denotes the inverse of the 3x3 matrix q'^S^^ + . 

In the non-relativistic hmit, is suppressed by powers oi v/c compared to II'^^ and this 
condition becomes simply g^ — n°°(0, q) < 0, which will happen for some q if — n°''(0, q) < 0. 
[A simple mnemonic for this condition is to recall that, for isotropic distributions, — n°'^(0, q) 
gives the squared Debye mass m^. So — n°°(0,g) < is analogous to having the 
wrong sign, which one might guess could indicate instability] This electrostatic instability 
is routinely discussed in plasma physics texts^ and is typically called the two-stream or 
Buneman [28] instability. We will refer to instabilities discovered by Condition 2 as "electric" 
instabilities, even though the situation is a little more complicated in the relativistic case. 

Because g*n°*(0, q) = 0, the q polarization does not contribute in the last term of (2.15). 
One may then restrict the sums over i and j to directions transverse to q and similarly 
project n to a 2 X 2 matrix in the transverse subspace. 

Were it not for the ie prescription in the retarded self-energy (2.8), n°*(0, q) would vanish 
by parity {p — >• —p). Using {x — ie)~^ = P.P. + iTi6{x), where P.P. denotes principal 
part, parity then gives 

n°^(0, q) = ine' [ v'5{v ■ q) q ■ Vpf{p), (2.16) 
Jp 

which is pure imaginary. We may therefore rewrite Condition 2 as 

q'-U'%0,q)-[lmW\0,q)] [g^ + n(0, q)] [Imn^°(0, q)] < . (2.17) 

If there are no magnetic instabilities, so that g^ + 11(0, q) is positive definite, then the last 
term contributes negatively to the left-hand side. In this case, a weaker sufficient condition 
for an electric instability is then that q^ — n°°(0, q) < 0. 

Our categorization of instabilities as "magnetic" or "electric" is not always physically 
significant. We use this terminology to tell us whether the existence of the instability was 
indicated by Condition 1 or Condition 2 respectively. Condition 1 is related to the transverse 
eigen- values of the zero- frequency self-energy n*-'(0,g), and Condition 2 is related to the q 
polarization [see Appendix A 2]. However, the actual growing unstable solutions to the 
dispersion relation (2.11) have non-zero frequency cu (with Imo; > 0). The corresponding 
polarizations can look quite different from those of 11(0, q) in generic situations where eigen- 
polarizations of W^{uj,q) are not fixed by symmetry. In particular, we will see a specific 
example in Sec. Ill where, for certain limiting cases of q, the actual growing unstable modes 
for a "magnetic" instability and "electric" instability end up being the same up to a rotation 
about an axis of symmetry. In the qualitative discussion that follows, however, we will focus 
on the Penrose conditions themselves and the nature of the self-energy at zero frequency. 



For example, a small sampling of references we have found useful is Refs. [25, 26, 31]. 
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C. Qualitative Origin of Instabilities 



We will now give a brief qualitative review of the origin of the instabilities we have 
discussed, in part because we do not know of a really good qualitative review in the non- 
relativistic plasma literature (especially of how the instability is saturated in the Abelian 
case), and in part to emphasize in a few cases differences for ultra-relativistic plasmas. In 
the context of the QCD plasma literature, a brief sketch is given by Mrowczyhski in Ref. 
[4], which we will expand upon. Understanding how the instabilities are generated is useful 
for intuition in understanding what types of instabilities may arise for a given f{p). 

We will start with a formal observation about Condition 1 for Weibel instabilities. It is 
useful to integrate the formula (2.8) by parts to obtain^° 



—CO + V ■ q — ie {—uj + v ■ q — ie^ 



+ f"'" (2.18) 

dp [—00 + vq — ley 



The last term vanishes for ultra-relativistic theories, for which v = 1 and dv/dp = 0. The 
interesting property of this expression is that, if one were to ignore the ie prescription, the 
integrand above would be positive semi-definite for a; = 0. If not for the contribution near 
the singularity, n*-'(0,q) could never have a negative eigenvalue. To see this, consider any 
unit polarization vector e and use (2.18) to write 



V 



[e'^(v ■ qf - 2(e ■ q){e ■ v){v ■ q) + q^{£ ■ vf\ 



ff^(0, q)e^ = e' - f{p) ^^^^^ ^ !\: ^' '\ (2.19) 



p [v ■ q — iey 

whose integrand, if not for the ie prescription, would always be non-negative by the vector 
inequality 

a^(6 • cf - 2(a ■b){b- c)(c • a) + b\a ■ cf > 0. (2.20) 

The conclusion is that, in the pcrturbative limit, the origin of the Weibel instability 
must come completely from hard particles with v perpendicular to q, which is when the ie 
prescription plays a role. All other hard particles are stabilizing. As we shall review, the 
destabilizing particles with v perpendicular to q are just the perturbative limit of particles 
whose motions in the q direction are trapped by the fields. The importance of trapped 
particles to discussions of instabilities, and their relation to singularities of integrals in 
perturbative calculations, has a long history starting with the somewhat analogous case of 
electric instabilities in one dimension [29, 30]. Note also that the effect of the stabilizing 
particles in (2.19) generically diverges as cos^ — > like d{cos9)/ cos^ 9, where 9 is the angle 
between p and q. 

To understand all this qualitatively, consider first a gas of non-interacting hard particles 
in zero field, as shown in Fig. 3a. By our parity assumption, there is no current j{x) 
associated with these particles because, for all the particles going in one direction, there are 



Here and henceforth, we impUcitly assume free hard particle dispersion relations are isotropic so that the 
direction of v is the same as that of p. Kinetic theory is sometimes used to describe the effects of hard 
modes in classical thermal field theories on a discrete lattice, where this assumption would not be true. 
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(a) (b) 

FIG. 3: (a) Parity opposites in a gas of free-streaming particles, (b) Small deviations from straight- 
line trajectories of untrapped (stabilizing) particles in the background of a small magnetic field 
B = BBySm^qz). The arrows show the direction of motion of positively charged particles, or the 
direction opposite to that of negatively charged particles. The direction (±?/) of the magnetic field 
is shown at the top, and the vertical dashed lines correspond to B{z) = 0. The direction (=bx) of 
the average current j{z) caused by the particles is shown at the bottom. 

just as many particles going in the other direction. Now turn on a small magnetic field B 
with wave- vector qr, which we shall take to point in the z direction. Take B to be in the ±|/ 
direction with B = Bsy sin(gz), as depicted in Fig. 3b, where By is the unit vector in the y 
direction. We'll take the vector field A to be in the ±x direction, with A = Ae^; cos{qz) and 
B — 'V X A. Magnetic forces from the small B field will make charged particles shghtly 
wiggle around straight-line trajectories. The small wiggles in direction will cause the x- 
component of the current to be larger in some places and smaller in others, compared 
to the value obtained from the straight-line trajectory. This is shown at the bottom of the 
figure. These currents set up a magnetic field that opposes the original magnetic field and 
so stabilize against its growth. 

The radius of curvature of particle trajectories in the plane perpendicular to the magnetic 
field is 

where pj_ is the momentum in that plane. It's easy to estimate^^ that Av^ ~ eAv/p so that 
the efi^ect described gives a contribution to of order e'^A J^fv/p. This contribution to 
the current is finite in the cos 6* — > limit and does not explain the small cos 6* divergence 
discussed earlier. 

Besides this variation in the magnitude of for each particle at a given position, there 



The field has roughly the same order of magnitude over a range of z close to half a period A/2 = n/q of 
the magnetic field. Simple geometry gives that the resulting small change in Vx over half a period is then 
of order vX/R, which in turn is ~ evB/pq ~ evA/p. 
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FIG. 4: The same situation as in Fig. 3b but for particles with small enough Vz that their z motion 
is trapped. 

is a second effect: the amount of time the particle spends at a given z coordinate changes, 
because the component Vz of its velocity also wiggles. The particles will spend more of their 
time in regions of z where Vz is smaller (which are the regions where is bigger). The 
slowing/speeding effect of the z motion therefore generates a contribution to the average 
current j{z) which is also in the directions shown in Fig. 3b. One may parametrically 
estimate that the variation Av^ of Vz has relative size of order^^ [^Vz)/vz ~ eAvv^/ {pvD, 
giving 



This gives the d{cos6) / cos^ 6 divergence, since Vz — vcosO in our choice of coordinates. 
This effect is also stabilizing. 

So far we have focused on the contribution from particles that stabilize magnetic fluc- 
tuations. Now consider particles whose initial velocities are very close to being orthogonal 
to q. The curvature of trajectories caused by the magnetic field will trap the z motion of 
such particles, as depicted in Fig. 4. The two trajectories shown on the left of the figure 
generate currents opposite those shown in Fig. 3, which therefore generate a magnetic field 
that adds to the existing one and contributes to instability. The contribution of the tra- 
jectory on the far right can be appreciated as follows. Imagine initially two particles, one 
at the point marked a and one at the point marked 6, have parity-opposite momentum in 
the ±a; directions. Then imagine suddenly turning on the magnetic forces in the problem. 
The upward moving particle at a will follow the middle trajectory and continue to give an 
upward contribution to j in roughly the same region of z. The downward particle at 6, 
whose contribution to j was originally canceling, will drift far away and so spend only a 
portion of its time close to its original z. That leads to the generation of a new upward 
average current in the region of their original z. 



Motion in a magnetic field conserves \v\. This means that Vx^Vx + Vz^Vz ~ 0, so that {Avz)lvz ^ 
{I^Vx)vx/v1. Then, use the previous estimate of Au^. 




(2.22) 
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(a) (b) (c) 

FIG. 5: Some qualitatively different examples of anisotropic distributions f{p) together with some 
specific examples of directions for q. The plots show the shape of surfaces of constant f in p space. 



For very small magnetic fields, only a perturbatively small fraction of phase space, cor- 
responding to very small If^l, is trapped. However, each such particle contributes to the 
current separation an 0(e) amount, unsuppressedhj the size of A. In contrast, the untrapped 
particles, making up almost all of phase space, each contribute a stabilizing influence; but 
the size of that stabilizing influence is perturbatively small (that is, suppressed by a power 
of the size of A). It is for this reason that trapped and untrapped contributions can be the 
same order of magnitude even in the perturbative limit. 

In the case of isotropic hard distributions f{p), the contributions of trapped (desta- 
bilizing) and untrapped (stabilizing) particles must cancel to give the isotropic result 
n*-' (0, q) = 0. This delicate cancellation allows us to understand the presence or absence of 
magnetic instabilities in certain simple situations. If we start with an isotropic distribution 
and, for a given q, remove some of the untrapped particles, then we must have magnetic 
instability. This is the case for the q shown for an oblate distribution in Fig. 5a. If we in- 
stead add untrapped particles, we must have stability, such as for the q shown for a prolate 
distribution in Fig. 5b. In 5c, assume that the number of particles outside the surfaces is 
insignificant. Then this is another example where there is no magnetic instability (for the 
q shown) because the number of trapped particles is insignificant. 

Finally, we will review the related picture for electric instabilities. We will focus on 
n'^°(0, q), which gives the charge response to a sinusoidal electric potential. Consider a small 
static electric potential A^ = (j)cos{qz) and the corresponding electric field E = — VA^ = 
Scz sin{qz). First think about the one dimensional problem — that is, particles moving in the 



Details depend on the distribution f{p) and the direction q. However, for generic distributions and 
generic directions q, one can parametrically estimate that trapping typically occurs in weak fields when 
R\cos0\ < 1. From (2.21), that is | cos 6*1 < y/eA/p, which indeed is formally arbitrarily small in the 
perturbative limit e — *■ 0. In a hand-waving sense, the ie prescription in the denominator of (2.19) can 
then be thought of as being of order \JeAfp. In any case, one can use (2.22) and integrate down to this 
I cos 1 to parametrically estimate the stabilizing untrapped contribution to be of order e fv / ^JeAjp. 
One finds a similar estimate for the order of magnitude of the trapped contribution, and the leading 
0(e^/^) pieces of these contributions cancel to leave the perturbative O(e^) result for n*-' (0, q). 
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FIG. 6: Qualitative picture of particle motion in a sinusoidal electric potential, (a) An untrapped 
non-relativistic particle in one dimension, (b) Untrapped and trapped particles in two or three 
dimensions. 

z direction, depicted in Fig. 6. If the field is small, typical particles will be able to move over 
the electric potential barriers. By energy conservation, however, positively charged particles 
will move slower at the maxima of and faster at the minima. So they will spend more 
of their time near the maxima, and so their contribution to the average charge density will 
be greatest there, which will generate a yet larger ^4° there. (Negatively charged particles 
will spend more time at the minima of A^, which will also enhance the magnitude of A^.) 
Untrapped particles are therefore destabilizing, which is the opposite of what happened in 
the magnetic case. In contrast, positively charged particles which have tiny velocities are 
trapped in the minima of AP and so contribute positively to the charge density there, which 
reduces ^4°. The electrically trapped particles are stabihzing. The one-dimensional situation 
is often discussed in the non-rclativistic case but is not relevant to ultra-relativistic particles 
because such particles always move with v = 1. However, in the three-dimensional case, 
can become smaller or larger as particles move through minima and maxima even while v 
remains the same, just as in our discussion of the magnetic case. Examples of electrically 
trapped and untrapped particles are given in Fig. 6b. 

Isotropic distributions are electrically stable and associated with a positive value of 
— n°°(0, q) given by the squared Debye mass m^. Small deviations from isotropy will not 
satisfy Condition 2 of Sec. II B, but significant deviations from isotropy can. 

An example where one can simply conclude that Condition 2 for electric instability is 
not satisfied is an oblate distribution, such as shown in Fig. 5a, with q pointing along the 
axis of axial symmetry, which we will take to be the z axis. For reasons similar to those 
given previously for the magnetic case, except that now the role of trapped and untrapped 
particles is reversed, — n°'^(0, e^) > 0. Axial symmetry and the Ward identity (2.9) imply 
n°'(0, Sz) vanishes, and so Condition 2 will not be satisfied in this situation. 
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D. Saturation of Instabilities 



When discussing the hnear response of the gauge fields, there is no significant difference 
between Abehan and non-Abehan theories. When the instabihties grow large enough to 
become non-linear, there is. It is useful to understand qualitatively how instabilities saturate 
in the Abelian case, in order to contrast it with the non-Abelian case. 

The saturation of Abelian electric instabilities is straightforward to understand qualita- 
tively. Return to Fig. 6 and imagine making the electric field arbitrarily large (holding q 
fixed) . Then essentially all particles will become trapped in the q direction near the bottom 
of the potential wells, which means there will be no instability. So, if there was originally an 
instability, the field will eventually grow large enough for it to saturate. This will happen 
very roughly when the electric potential energy eA^ becomes of order the typical energy 
associated with the motion of the particles in the q direction. In the non-relativistic limit, 
that is 

pfi pIq. 

~ , ^-47' (2-23a) 

where M is the particle mass and p\\ indicates the typical value of the component of p parallel 
to q. In the ultra-relativistic limit, 

2 2 

^0 ^ ^ , £; ~ M . (2.23b) 
ep ep 

This saturation was used long ago to construct non-linear wave solutions in the non- 
relativistic Abelian case, known as Bernstein- Greenc-Kruskal waves [30]. 

To qualitatively understand the saturation of magnetic instabilities, return to the pictures 
of Figs. 3 and 4, but now make the magnetic fields arbitrarily large (holding q fixed). In 
this case, the typical radius of curvature (2.21) will be small compared to the distance in z 
over which B{z) changes sign. Instead of looking like Fig. 3b or 4, typical trajectories will 
instead look like Fig. 7. The particles follow small, nearly circular orbits in the magnetic 
field which slowly drift in the direction indicated due to slightly higher curvature on the side 
where B is slightly larger. This sets up currents, shown at the bottom of the figure, which 
create magnetic fields that oppose the original field. Essentially all particles are therefore 
stabilizing in the high field limit. So, once again, if there was originally an instability, the 
fields will eventually grow large enough for it to saturate. That will happen very roughly 
when the radius of curvature (2.21) becomes of order the wavelength A ~ 1/? of the field. 
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FIG. 7: Trajectories in an arbitrarily large magnetic field. The size of the orbits and the rate of 
their drift has been exaggerated. 



B^?ll, A^^. (2.24) 
e e 

A simple mnemonic for cases where the hard particles can be described as excitations in 
a quantum field theory is to realize that, when acting on hard excitations of momentum 
p, the gauge term in the covariant derivative D = d — ieA cannot be treated as a small 
perturbation to the derivative term when (roughly speaking) A> p/e. 

In contrast, if different unstable soft field modes with momenta of order q can interact 
with each other, then linearity breaks down when 

Ar^^ (2.25) 
9 

instead of ^4 ~ p/g, where we have switched notation from e to as a way of emphasizing 
that this is special to the non-Abelian case. We will see a specific example in Sec. IV. 



This and (2.23) represent estimates in generic cases of significant anisotropy. If, for example, the anisotropy 
is very tiny or q is very close to its maximum unstable value, then a very tiny higher-order correction to the 
perturbative calculation might change instability to stability. A brief and somewhat different discussion 
of rough estimates of the saturation condition in the non-relativistic case for both magnetic and electric 
instabilities may be found in Ref. [33], as well as references therein. There the saturation condition is 
estimated in terms of time rather than distance scales by equating the "bounce" time for an oscillation 
of a trapped particle to the linearized growth rate of the instability. For generic situations and generic 
unstable q, this is roughly equivalent to our rough conditions (2.23) and (2.24) if one assumes that the 
growth rate 7 is of order vq (where we will no longer try to distinguish between and v±). See, for 
example, Eqs. (1), with the condition wb ^ 7, and (114), with k small, of Ref. [33]. The estimate vq 
(which is <IC 9 in the non-relativistic limit) can be roughly understood by thinking of increasing u = ij 
from zero and noting that it first has an effect in the self-energy when w ^ vq, so that it significantly 
affects the integration over the denominators in (2.18). 
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Even for the Abelian case, the picture of Fig. 7 is a bit simplistic. It assumes that only 
one unstable mode with a given q has been excited. In most situations, an entire spectrum 
of unstable modes will grow together and then interact with each other once they become 
non-perturbatively large according to (2.24). The ideahzed case of a single Abelian mode has 
been studied in the non-relativistic literature and certain non-linear wave solutions have been 
found, known as magnetic Bernstein- Greene-Kruskal (BGK) waves [32, 33]. The generation 
of these non- linear waves has also been investigated numerically [33-35], in some cases with 
only a single unstable mode significantly excited and in other cases with many. 

III. INSTABILITY OF THE PLANAR MOMENTUM DISTRIBUTION 

We will now focus in detail on the planar momentum distribution of (1.1): 

fip,x) = Fip^)6ip,). (3.1) 



A. Magnetic instabilities 



UJ 



1. u} — magnetic stability analysis 

To check for evidence of magnetic instabilities, first consider the self-energy (2.18) at 

pip 



0, 



m^(o,4) = e^ 



p - q — ie (p - Q — i^)" 



(3.2) 



The planar distribution f{p) has support only in the PxPy plane and is independent of the 
direction of p in that plane. We may therefore factor out the angular dependence of the 
integral and write 



m^(o,g) 



pq-ie (p-g-ie)Vpeplane 



(3.3) 



Because of the axi-symmetry of the problem, we can assume without loss of generality that 
q points in the xz plane. Then Sy is an eigenvector of 11'-' because of y — > —y reflection 
symmetry. For a; = 0, the direction q is also an eigenvector of 11'-', corresponding to zero 

eigenvalue. Let n be the orthogonal direction in the xz plane, as depicted in Fig. 8. Let 9 
be the angle q makes with the z axis and (f) the angle p makes with the x axis. Then, for 
sin^^ 7^ 0, the eigenvalues of n*-' (0, q) are easily computed from (3.3) to be 



n^^(o, q) 
n""(o, q) 



0, 



27r 



27r 



n^2/(0,q) = m; 



r^" d<j) 

Jo 2^ 



1 + + 



1 + + 



cos^ 6 



sin'^ 



(3.4a) 
(3.4b) 



sm 



sin 9{cos(p — iey 



sin' 9 



^-mi^ cot' 9. 



(3.4c) 
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FIG. 8: Convention for labeling directions with respect to g, which is taken to lie in the xz plane. 



Only W^iQ^q) is negative.^^ By Condition 1 of Sec. II B, there is a magnetic instability 
associated with A polarized in the y direction (and so B in the h direction) whenever 

q < gmax(6') = rrioo cot 9. (3.5) 

It is easy to understand qualitatively why there is no magnetic instability associated 
with the n polarization of A, which would correspond to the magnetic field B in the ±y 
direction. Magnetic instability is caused by trapped particles, which perturbatively arc those 
with p orthogonal to q. For the planar distribution with q in the xz plane, the only trapped 
particles are those with p pointing in the ±y direction (for sin^ 7^ 0). If B is also in the 
±y direction, then the motion of these particles is not affected by the magnetic field, and 
so there is no destabilizing contribution. 



2. Magnetic instability growth rate in the ultra-relativistic limit 



We can get simple expressions for W^{uj,q) for general to if we speciahze to the ultra- 
relativistic limit, which wc will now do. In this case, it is straightforward to keep the 
frequency dependence in the manipulations that led to from (2.18) to (3.3) for the self- 
energy, giving 



ff^(a;,q) 



5'^ - 



-r]+p - q-ie 



+ 



{-rf + l)py 
(-77 + p-q-ie)Vpg plane 



where 



V = - 

1 



(3.6) 



(3.7) 



Since there is no suggestion of magnetic instabilities for other polarizations for sin 6* 7^ 0, let 
us focus on n^^. The nice thing about the simple planar distribution is that one can find a 
simple closed form solution from (3.6). Performing the integrals, one finds 



n^V,qr) = m^a + 



sin'^ 



sin^^ 



-1/2 



(3.8a) 



This instability corresponds to poles in Aa of Ref. [17] in the ^ +00 limit. The possibility of instability 
in this mode was overlooked in Ref. [6] . 
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FIG. 9: The growth rate vs. q for magnetic instabihties for the planar distribution for various 
values of 6. Prom the top down, the values are sin^ = 0.01, 0.1, 0.3, 0.65. The small dashed line 
at the top of (b) shows the sin — limit. 



The corresponding dispersion relation is 



(3.8b) 



By re-arranging terms algebraically and squaring, one can convert this into a cubic equation 
in ?7^, a subset of whose roots are the solutions to (3.8b). However, the explicit closed form 
of that solution is no more enlightening than a numerical solution. In Fig. 9a, we show the 
result for the growth rate 7 = Ima; = gimr/ of the unstable solution as a function of q for 
various values of ^. It is clear that the largest growth rates correspond to sin^ — > 0, which 
is not surprising since, for ^ = 0, all particles in the planar distribution become trapped 
particles, contributing to destabilization. Plotting 7/moo vs. q/qmax in Fig- 9b, with gmax 
taken from (3.5), it is clear that there is a limiting behavior as sin^ 0. For small sin^, 
the dispersion relation (3.8b) can be more compactly solved to give 



7 



rrir. 



fi 



rrir. 



+ /2 



1 

71' 



where 



h{z) 

f2{z) 



1/2 



(3.9) 

(3.10) 
(3.11) 



(2 - ^2)1/2 ■ 

Note that moo — ?maxSin^ for small sin^. fi{z) runs from to I/V2 and f2{z) from 1/^2 



to as 2; runs from to 00. For q <^ moo, one gets 7 ~ g in this small sin^ limit. The 
maximum growth rate for a given small sin 6* occurs at ~ ^71.^(1 sin^ Oy^^. 
For sin 6* = 0, the result becomes simply 



7 



fi 



rrir 



(sin^ = 0). 



(3.12) 



which, as g — s> cxD, approaches the maximum growth rate for the planar problem of 7 = 
moo/ \/2- By symmetry the behavior of fl"^^ must be the same as that of 11^^ at sin ^ = 0, 
and so there are two rather than one magnetically unstable modes at exactly sin ^ = 0. 
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FIG. 10: Qualitative behavior of n"'"'(0,1g') and 11^^(0, q) vs. q for a planar distribution with small 
thickness A^. The curves for sin^ 3> depict the formulas of (3.4). 



3. Thick planar distributions and related matters 

One can now ask what happens if the planar distribution has a tiny thickness, so that 
S{pz) is replaced by something with a small but finite width Ap^. Let ~ {Apz)/p be 
the angular width of the distribution, where p is the momentum scale that dominates /(p), 
and assume A^ <^ 1. The effect this will have is to smear out the 9 dependence of previous 
results over A^^. In particular, 11^^(0, q) will now be finite for q in the ±2; direction, with 
(3.5) replaced at ^ = by 

qmU0=O) - ^ ■ (3.13) 

Also, there will now be two unstable modes for sin^^ < A^^ rather than just at sin^ = 0. 
The qualitative behavior of II'^'^(0, q) and 11^'^ (0, q) as a function of 6 is shown in Fig. 10. 

Though the preceding is a correct analysis of "magnetic" instabilities as defined by Con- 
dition 1 of Sec. II B, it is qualitatively a bit misleading. Instead of focusing on the zero 
frequency self-energy used in Condition 1, let us make some very rough qualitative argu- 
ments about the self-energy at the non-zero frequency u = i'-f corresponding to the growing 
unstable mode. The origin of the formal divergence of the result (3.4c) for 11^'^ (0, q) as ^ 
is in the denominators {—f] + p- q — ie)^ for the zero-thickness (A^ = 0) planar distribution. 
For rj = 0, this is negative (destabilizing) only when p ■ q = 0. For r] = i'j /q, it is instead 
destabilizing when p ■ q ^ 'y /q- For q close to the z axis, the average value of {p ■ q)^ is 
of order The qualitative effect of non-zero 7 is therefore similar to thickening out the 
planar distribution by A9 ^ j/q. For small 9, we have found magnetic instabilities with 
7/g ~ 1 when q <^ ^max, which is therefore similar to A^ ~ 1. We therefore expect that for 
small ^ ^ 1 and q <S ?max, the smearing due to non-zero 77 should be enough to allow an 
x-polarizcd unstable solution that looks just like the y-polarized solution we have already 
found. Indeed, we shall discover just such a solution in the next section, where we turn to 
the "electric" instabilities associated with Condition 2 of Sec. II B. 
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B. "Electric" instabilities in the ultra- relativistic limit 



1. ui = electric stability analysis 



Let us look for electric instabilities by checking Condition 2 of Sec. II B. First consider 
n°°(0, q). Starting from the general expression (2.8) for the self-energy and integrating by 
parts, one finds 



n0°(0,q) = e2 / f{p) 
Jp 



vp 



{p q- ief 



^ 1 dv 
dp 



(3.14) 



Let's focus on the case sin 9 0. For the planar distribution function, the angular integral 
of {p • q — ie)~^ over p in the xy plane gives zero, and one finds the q-independent result 



n0°(0,q) = e2 / f{p) 

J n 



■ 1 I dv 
vp dp 



(3.15) 



This result gives zero in the non-relativistic limit. We will focus on the ultra-relativistic 
limit, where 

n°°(0,q) = -m^, (3.16) 

which is negative and so stabilizing. 

But now we need to examine n°*(0, q). From y — > —y reflection symmetry, one has 
n°^(0, q) = and from the Ward identity, n°*^(0, q) = 0. It remains to calculate n°"(0, q). 
One can use (2.16), but we find it simpler to proceed similarly to earlier calculations with 
the planar distribution. Integrating (2.8) for the self-energy by parts. 



U'\0,q) = e' 



f{p) -q'v-q 
p p (v-q- ie)2 



(3.17) 



In particular, we can write 



n""(0, q) = e' 



fip) 



'— cos 9 cos I 



i cos 9 



p \{sm9 cos4> — ie)"^ / ^ sin^ 6* 
In the ultra-relativistic limit, this is just 

cos^ 



p 



(3.18) 



n""(0, q) = -im, 



°°sm^9' 



(3.19) 



Now we can apply Condition 2 of Sec. II B, which indicates an electric instability if 

- n™(0, q) + nO^(0, q) [q' + U(0, q)]'' W%0, q) < 0. (3.20) 



Making use of the result (3.4b) for n"'"'(0, g), this condition becomes, for the ultra-relativistic 
planar distribution. 



2 2 ml,cot^6' 



q^ sin 9 -\- 



(3.21) 
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The left-hand side is minimized for g = 0, in which case the inequahty gives 1 — cot^ ^ < 0. 
So, there will be some mode of electric instability whenever q lies within 45 degrees of the 
z axis.^^ For a given 9, the maximum q which gives an instability is given by (3.21) as 



cos e (4 + 008^^)1/2 - 



sin^ 9 - 1 



.1/2 



5max 



'OO 



2 sin^ 9 



(3.22) 



Since — n°°(0, q) > 0, which is stabihzing, the "electric" instabihty in this case is not due 
to the simple picture of electric instabilities discussed in Sec. II C but instead depends on the 
coupling of charge and current fluctuations through 11*°. An example of a distribution with 
an electric instability whose origin does have the simple interpretation of Sec. II C [having 
n*°(0, q) — 0] is given in Appendix C. 

2. Electric instability growth rate 

In Appendix D, we briefly summarize the calculation of the dispersion relation for polar- 
izations of A in the xz plane {qn plane), just as (3.8) gave the dispersion relation for the y 
polarization. The result is^^ 



where we have introduced the dimensionless variable Q = qlm^o- The unstable solutions are 
plotted in Fig. 11 for various choices of 9. Unhke the magnetic case, there are no solutions 
for angles more than 45° from the z axis. For any given the maximum instability occurs 
for sin^ 0, which is a limit in which the unstable solution to (3.23) becomes 



which is just the same as (3.12) for the magnetic case. One may check (see Appendix D) 
that the corresponding polarization of E is in the x direction. This is just the solution 
previously predicted for small sin^ in Sec. Ill A 3. 

One may also explore taking sin 6* — > for g ~ gmax as in Fig. lib (rather than q fixed). 
In this limit, both ^max and the functional form for 7 in terms of q/qxaas. are different than 
in the magnetic case. We do not have a simple closed-form expression for this limit in the 
electric case. 

Summarizing our results for all types of instability (electric and magnetic) for the planar 
distribution (1.1), we have found that there are two polarizations of instability for q within 
45 degrees of the z axis and one polarization otherwise. The most unstable modes are those 
with q close to or along the z axis. 



This instability corresponds to poles in of Rcf. [17] in the ^ +00 limit. 

By re-arranging terms algebraically and squaring, one can convert this into a quintic equation in 77^, a 
subset of whose roots are the solutions to (3.23). 





1/2 



l-g^sin^^ =0, (3.23) 




(3.24) 



23 



2 4 6 8 10 

(a) 



0.2 0.4 0.6 0.8 1 

(b) 



FIG. 11: The electric instability growth rate vs. g for the planar distribution for various values of 
d. Prom the top down, the values are sin0 = 0.01, 0.1, 0.3, 0.65. The small dashed line at the top 
of (b) shows the sin^ — limit. There are no solutions for sin^ > l/\/2. 

IV. NON-ABELIAN EXAMPLE OF SATURATING THE INSTABILITY 

We now give a concrete example of one particular way in which non-Abelian interactions 
between soft fields can halt the continuing growth of the amplitude of an unstable mode. 
Specifically, we will show that there exist static non-linear wave solutions to the soft field 
equations of motion which, in the linear approximation, would correspond to unstable modes. 
The soft field amplitudes will be of order qjg^ as in (2.25). We will assume that the soft 
modes of interest have momentum q much smaller than the momentum p of typical hard 
particles, so that gAj]) <S 1. This means that the interactions of soft fields with each other 
must be treated non-perturbatively, but the interaction of soft fields with individual hard 
particles is perturbative. Wc can therefore ignore non-linear interactions induced by the 
hard particles and focus just on the perturbative self-energy 11 they induce for the soft 
fields. The relevant soft field equations are then the non-Abelian Maxwell equations with 

given as before by the perturbative self-energy: 



where F^^ is the full non-Abelian field strength and D is the covariant derivative d — igA 
in the adjoint representation. The only assumption we will make about the hard particle 
distribution is that it is axi-symmetric about the z axis so that (4.1) also has this symmetry. 

A. A solution 

Consider any SU(2) subgroup of the gauge group. Within this subgroup, we take the 
group structure constants to be the SU(2) ones /"^"^ = e"'"^. Working in Lorentz gauge 
{df^A'^ — 0), our non-Abelian Maxwell equation is then 



D^F^"' ~ WA, 



(4.1) 



nA^"" + ge' 



a 



(4.2) 



We will look for solutions with the ansatz 




ARaniqz), a = 1 or 2, and = 1 or 2; 
otherwise. 



(4.3) 
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where a is the adjoint color index and Rij{0) is the 2x2 rotation matrix 

n f cos6' sin6'\ 

^-(^)=(,-sin^ cos^j- (4-4) 

The derivative in the second term of (4.2a) then vanishes. Using RuRdi = one then 
obtains 

{-q^ - 2g^A^)ARaAqz) = n^(0, qe,)ARauiqz) , (4.5) 
where = = IP^. We have a solution if 

A'^-^[q' + U40,qe,)] (4.6) 

is positive. That will occur exactly when + 11^(0, ge^) < 0, which is precisely Condition 1 
of Sec. 11 B for a magnetic instability associated with small-amplitude waves with the same 
q = qSz- Note that in our static non-linear solution, A oc 1/g and so is non-perturbatively 
large. 

The energy density stored in the soft field is ~ q'^A^. Since A ^ q/g and q ~ moo, the 
energy density in the soft field is ~ m^/g'^. This is the same energy density we would obtain 
if all potentially unstable modes (of order half of all modes with q < m^) carried occupation 
number ~ l/s'^- The latter is the situation we actually expect, when the instability saturates. 



B. A simple analog in scalar theory 

To motivate the ansatz (4.3), it is instructive to note that there are analogous static 
solutions of simple 0^ theory of a complex scalar (j) with a Mexican hat potential 

y(0) = VI0r + |A|0r- (4.7) 

Here the equation of motion is 

□0-A|0|^0 = -//^0, (4.8) 
and — plays the role of the destabilizing negative 11. By taking the ansatz 

(t>{x) = Pe^«^ (4.9) 
where e*^ is hke the rotation matrix R{6), one obtains 

V'^-\[q'-l,\ (4.10) 

which gives a static solution for q < fi. These g's correspond to what would be the unstable 
modes in a small fluctuation analysis about = 0. 

The simphfications attending the use of a "helical" ansatz that rotates with 9 — qz, 
like the two ansatzes above, were also used in Ref. [32] to find a special case of non-linear 
solutions in the (qualitatively different) Abelian gauge theory problem. 
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As far as we know, the scalar solution, like the non-Abelian one preceding it, is not 
relevant to any realistic physical situations. In typical applications where unstable fluctu- 
ations are growing, one will have an entire spectrum of different growing modes which will 
interact with each other non-linearly once they grow large enough. In the scalar problem, 
we would not expect a random spectrum of initial small fluctuations about (6 = to grow 
into a solution like (4.9). Instead, we expect this system to equilibrate, dissipating its energy 
into small random fluctuations about a minimum of the potential energy. Similarly, though 
our non-Abelian solution gives an example of how non-Abelian interactions can saturate the 
growth of instabilities, we expect the actual situation to be much more complicated. 

As a flnal note on how things can be very different when only a limited number of modes 
are involved, consider a different ansatz we could have considered for the gauge theory 
problem that only considers modes with a single (adjoint) color: 

A;ix) ^ A{x) 5^%. (4.11) 

With this ansatz, the non-Abehan non-linear equations (4.3) reduce to the Abehan linear 
equations (2.7). In this unnaturally idealized situation, the growth of instabilities would 

not be halted by soft interactions but, just like the Abclian case, would continue to grow 
until there were non-perturbative effects on the hard particles. We suspect that such a 
conflguration in a non-Abelian setting would be unstable to the formation of large flelds in 
other color directions, but we have not shown this explicitly. 

V. CONCLUSION 

We have seen that QCD plasma instabilities play an important role in the idealized 
theoretical limit of high energies and small coupling constants and drastically modify the 
bottom-up thermalization scenario. We have also seen reasons why the saturation of these 
instabilities should be qualitatively different from the Abelian case. An important goal for 
future work is to understand what replaces the original bottom-up scenario and determine 
in detail how these instabilities parametrically affect the rate of equilibration. 
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APPENDIX A: SUFFICIENT CONDITIONS FOR INSTABILITY 

In order to analyze the existence of instabilities, it will be convenient, initially, to phrase 
the discussion solely in terms of II*-'. Note that II^*^ and 11°' can be converted using the 



For q> n, one can alternatively find no7>static solutions of the scalar theory of the form (f){x, t) = "Pe"^*. 
As far as we know, these arc not useful for anything either. A broader class of related solutions for the 
quark-gluon plasmas has been investigated for stable, equilibrium plasmas in Ref. [36]. 
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the Ward identity (2.9). We will follow Mrowczynski and others by rewriting the linearized 
Maxwell equation (2.7) in the form^^ 

D'^iu, q) = [(-a;' + q^'^ - q'q^ + U'^{u;, q)]E^ = 0. (Al) 

One quick way to obtain this is to take the time derivative of (2.7) in = gauge with 
jjL — i. The question of the existence of instabihties becomes whether D^^{uj, q) ever has a 
zero eigenvalue for some uj with Imo; > 0. 



1. Condition 1 

In searching for solutions to the dispersion relation with complex cj, it will be convenient 
to drop the ie prescription in the self-energy (2.8) and simply absorb it into the value of ui: 



—UJ -\-v ■ q 



(A2) 



Changing integration variables p — > —p in this expression then implies W^{uj) = W^{—uj). 
Complex conjugation then yields^° 

m^(cu) = [rf^(-cu*)]*. (A3) 

In particular, consider pure positive imaginary uj, and write ou — i'j with 7 real. Then (A3) 
implies is a real matrix. Since it is also a symmetric matrix, it therefore has real 

eigenvalues. 

The expression (A2) for W^{iu, q) is bounded in magnitude as a; — > ioo. In this hmit, D^^ 
given by (Al) then becomes 

D'^ {u, q) -u^S'^ = -i^5'\ (A4) 

with eigenvalues becoming identical and going to +00 as 7 — > 00. On the other hand, for 
7 = 0, 

D^i(0, q) = q^5'^ - q\^ + n*^ (0, g), (A5) 

which is the 3x3 matrix of Condition 1' of Sec. II B. Suppose D^^{0,q) has one or more 
negative eigenvalues, such as shown in Fig. 12a. By continuity of the real eigenvalues of 
£)*^ (0, i^) as 7 is varied from to 00, there must then be some positive 7 for which D^^{0, i^) 
has a zero eigenvalue, which corresponds to an unstable solution. And there must exist one 
such solution for each negative eigenvalue of D'-'(0, g). This proves Condition 1', which is 
equivalent to Condition 1. 



^'^ This equation is gauge invariant in Abelian theory. Gauge invariance in non-Abehan theory is not manifest 
simply because we are studying instabilities by linearizing in the fields. The equation is gauge invariant 
up to non-linear terms that have been neglected. 

^° For w = ct^o + ie with t^o real, this implies W\^^^{ujq) — — n!^dv(~'^o) = ^[nret(~'^o)]*> where "ret" and 
"adv" denote the self-energies for retarded and advanced Green's functions. 
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(a) (b) 

FIG. 12: Examples of the continuity of eigenvalues of D^^{i^, q) with 7, depicting (a) one magnetic 
instability [Condition 1], and (b) one electric instability [Condition 2]. 

2. Condition 2 

There is always a zero eigenvalue of £)*-' (0, q) associated with the longitudinal polarization 
q, and the above argument does not tell us whether there is a 7 > solution D^^[i^, q) — 
associated with the continuous evolution of this eigenmode as 7 varies from to 00. However, 
suppose we knew that the corresponding eigenvalue was negative for very small positive 
7. The previous eigenvalue continuity argument would then guarantee a corresponding 
instability, as depicted in Fig. 12b. 

Let us therefore investigate the small eigenvalue A of D'^^{i'~f,q) for arbitrarily small and 
positive 7. Let T represent the two directions transverse to q. Making use of the Ward 
identity (2.9), the matrix D^^(uj, q) for small oj has the form 



D{uj,q) 



n'^u^q) q' + U-^^{u;,q)J \ |nT°(0,qr) q' + U^^{0,q)^ 

(A6) 

where 11(0, q) above should be understood (as in the main text) as the retarded self-energy 
at zero frequency, Il{ie, q). In the limit of small u, the small eigenvalue of such a matrix is 

A ^ ^ [-q' + n-^" - n°^(,^ + n^^)-^n°^_^^ , (A?) 

where [g^ + 11^^]"^ is the inverse of the 2x2 matrix + 11^^. Condition 2 of Sec. II B is 
just the condition that the eigenvalue (A7) be negative for small imaginary cu. 

3. Relation to Nyquist analysis 

In textbooks [25, 26], it is more usual to derive the Penrose criteria from a Nyquist analysis 
rather than to use continuity arguments as above. However, the usual analysis becomes more 
complicated than textbook examples for general distributions f{p) and general directions 
q, and we find the continuity argument simpler to implement for reasons we now explain. 

The textbook derivations work whenever one can use symmetry to reduce the matrix dis- 
persion relation (Al) to single-component dispersion relations. If one is looking for solutions 
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(a) (b) (c) 

FIG. 13: Examples of Nyquist contours D{C) in the complex plane. The large outer circle repre- 
sents the image of the semi-circle at infinity, and the rest is the image of ^ + ie. [Diagram (c) is 
not just a fanciful example. It is also the form of the Nyquist diagram associated with the function 
rf~'^ times the left-hand side of (3.23) for certain values of q and sin^ that give instability. The 
single winding denotes a single associated unstable mode.] 



to a single component relation of the form 

D{u;, q) = -up" + + n(a;, q) = 0, (A8) 

and if 11 is an analytic function of uj in the upper half plane, then one can count the number 
of zeros of D{iS) for fixed q (which gives the number of unstable solutions) using the Residue 
Theorem as 

Jc 2ni D du} J 27ii ' ^ ^ 

where the contour C encloses the upper half uj plane by (i) running infinitesimally above 
the real axis (which we will denote 3? -|- ie) and then (ii) closing in a semi-circle at infinity. 
The last form above shows that N is just the winding number about the origin of the image 
D{C) of the curve C. For |a;| — >• oo, the self-energy given by (2.18) is bounded, and so 
D{uj) — t^^, so that the semi-circle at infinity maps to a circle at infinity that almost 
closes. Fig. 13 shows various examples of possible images D{C). The winding number in all 
these cases is even if 11 (ie) is positive and odd otherwise, for the following reasons. Count 
the winding number by counting the number of times (and sense) that D{C) crosses the 
negative real axis. The curves D{C) are symmetric with respect to complex conjugation 
because the conjugation property (A3) of the self-energy implies 'n.{x + ie) = U* {—x + ie) for 
real x, so that the image of the half-line infinitesimally above the positive real axis must be 
the conjugate of that above the negative real axis. This conjugation property means that, 
except at UJ = ie, any crossing of the real axis by D{C) must be at a point where the real 
axis is crossed twice. As a result, all the crossings except for oo = ie do not affect whether 
the winding number is odd or even, which is solely determined by Il{ie). If Il{ie) is odd, one 
knows that the winding number cannot be zero, which proves the sufficiency of the Penrose 
criterion. 

The difficulty with the above analysis in the multi-component case is that it relies on D{uj) 
being analytic in the upper half plane. Suppose that in the general multi-component case 
of (Al) one wants to develop three independent instability conditions, corresponding to the 
three eigenvalues of the matrix D. One could imagine algebraically solving for the eigenvalues 
to get three one-component equations. However, solving the characteristic equation to obtain 
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eigenvalues introduces root singularities and branch cuts. If any of these appeared in the 
upper-half plane, it would destroy the analyticity assumption of the Nyquist analysis. It 
may be possible to construct an analysis of this sort, but we found it easier to avoid the 
issue altogether. 

One may forgo finding independent conditions and simply calculate the total number of 
unstable solutions by doing a Nyquist analysis with the single equation det D = 0. However, 
there are not useful corresponding Penrose criteria because the sign of det D{ie) would 
depend on whether there were an even or odd number of unstable modes. 

For cases where the problem can be reduced to one-component equations like (A8), one 
often sees the Penrose criteria used as necessary as well as sufficient conditions for instability. 
This happens, for example, in situations where one can argue that the sign of Imn(a; + ie) 
is the same as the sign of x for real a; 7^ 0, so that the only place DiJR. + ie) can cross the 
real axis is for cu — ie. Then the winding number is either exactly zero or one, as in Figs. 
13a and b, depending on the sign of Il{ie). For example, consider the case with q in the z 
direction, f{p) = F{pz)H{p±), and polarization in the x direction. Taking the imaginary 
part of (2.8) for real uj gives 

lmW%uj,q) = e'7r J ^ {v,f5{-i^ + v,) , (AlO) 

where 77 = oj/\q\. The sign of this result is then just the sign of dF/dpz evaluated at Vz — rj. 
If F{pz) is a monotonically decreasing function of \pz\, then the sign of the result will indeed 
be the sign of cu, and the Penrose condition would be necessary as well as sufficient for that 
mode. 

Unfortunately, matters are more complicated when q does not point in a symmetry 
direction. Suppose again that f{p) had the form F(pz)H(p±), with both F(pz) and H(p±) 
monotonically decreasing functions of \pz \ and p±. One could easily make a similar argument 
as above if Imn*''(u;,q) were a positive definite matrix for positive real u. One can check 
positive definiteness by checking e^lmW^e^ > for arbitrary real spatial polarizations £. 
Eq. (13) gives a result of the form 

£jlm[n'^'(a;,qf)]£j = / [positive] 6 {-rj + v • q) q -Vpf. (All) 
Jp 

Unfortunately, q ■ 'Vpf is not necessarily positive when v ■ q is positive, as shown in Fig. 
14. We are unaware of a general proof that the Penrose criteria are necessary conditions for 
instability in such cases. 



APPENDIX B: ANISOTROPY CREATES INSTABILITY (CONDITION 1-B) 



First, we will show that non- vanishing W^{0,q) implies magnetic instabilities. Start by 
considering the trace of the spatial part of (2.12), and then spatially average that trace over 
the direction of q: 



{W\0,q))^ = e' I 
Jp 



dfjp) 
dp^ 



V ■ q — ie 



(Bl) 
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q 



FIG. 14: An example (shown in the xz plane) of /(p), p and q for which v ■ V/ and v ■ q do not 
have the same sign. The ellipses indicate curves of constant /. 

It is easy to see that 



because the left-hand side (i) gives unity when dotted with v, and (ii) must be proportional 
to V by rotational invariance. We therefore obtain 



But the trace 11" is the sum of the eigenvalues of ll*-'. What we have shown is that the 
eigenvalues of n*-' (0, q), averaged over all eigenvalues and all directions of q, average to zero. 
Therefore, either all eigenvalues in all directions are identically zero, or some eigenvalue in 
some direction must be negative (implying an instability). In particular it is a necessary 
condition for stability that the trace 11" (0, q) should vanish for all q. 

Now we show that anisotropic M.{p), defined by (2.13), implies non-vanishing trace 
n"(0, q) for some q and therefore implies the existence of instabihties. Taking the trace and 
setting a; = in (2.18), we have 



Note that the factor in brackets docs not depend on the magnitudes p and f of p and v but 
only on their common direction p. Performing the integral over p gives a factor of M.{'p) 
defined by (2.13): 




(B2) 



(ff^(0,q)>. = 0. 



(B3) 




(B4) 




(B5) 



Now decompose M. into spherical harmonics: 




(B6) 



Im 
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By rotational invariance, Eq. (B5) will generate an angular dependence proportional to 
Yim{q) for each term proportional to Yim{p)- So 



n"(0,g) = ^KiaimYimiq), 



Im 



where 



1 + 



1 



Yimip) 



(B7) 



(B8) 



does not depend on m, again because of rotational invariance. 

The only way that (B5) can vanish for all q is if all of the aim^i vanish. If Ai{p) is 
anisotropic, then one of the aim must be non-zero for / > 0, and it must be an even / 
because of our universal assumption that f{p) is parity symmetric. If we knew that ni ^ 
for all even Z > 0, we would then know that 11** (0, q) cannot vanish for anisotropic M.{p). 
Therefore we turn to the evaluation of ki. 

To evaluate k,i, first choose m = in (B8). Next, use Wigner D functions to re-express 
the spherical harmonic Yio{p) defined with respect to a fixed z axis in terms of spherical 
harmonics y^^) (p) defined with respect to the direction q: 



YUp)-J2D'mm'iQ)y£kp)- 



(B9) 



Let 9 be the angle between p and q and note that 



1 + 



(cos^ — iey 



(BIO) 



vanishes unless m' = 0. So 

1 



1 + 



{pq- i^y 



Yioip)) =DUq) 



1 + 



{pq- i^f 



Yioiq) 



1 + 



Yifip) 
1 



(cos^ — ie)2 



P;(C0S^)> , (Bll) 



which gives 



1 



dx 



-1 



1 + 



{x — iey 



where Pi{x) is the l-th Legendre polynomial. This integral gives' 
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(B12) 



10 



III 



r(| + i) {i-m 



(B13) 



for even I, which indeed does not vanish for any even I > 0. 
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This may be derived by an appropriate analytic continuation of Eq. (7.126.1) of Ref. [37]. 
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This completes the proof that anisotropic A4{p) impUes magnetic instabilities. What 
about the converse? Suppose Ai{'p) is isotropic. In the ultra-relativistic limit (f = 1), we 
can perform the radial p integral of the expression (2.18) for the self-energy to write 



n 



M{p) 



5'^ 



_|_ 



+ 



{—if + l)p^fp 



with 77 



-r] + p-q-ie [-r] + P ■ Q ■ 
uj/q. II M.{p) is isotropic, then we can factor it out to give 



(B14) 



Wi{u,q)^ml,l5^^- 



-T] + p - q 



+ 



le 



-T] + p ■ q — iey 



(B15) 



All of the details of the distribution f{p) have factored out into the single normalization 
constant m^. Except for this normalization, there is no difference between the general case 
of isotropic A4{p) and the specific case of equilibrium distributions. In equilibrium, there 
is no instabihty, and so there cannot be one for any isotropic (One may check this 

explicitly using the standard results for the equilibrium self-energy [38, 39].) 



APPENDIX C: INSTABILITY OF THE ULTRA-RELATIVISTIC LINE MO- 
MENTUM DISTRIBUTION 

In the bulk of this paper, we have focused on instabilities associated with the planar 
(extreme oblate) distribution of (1.1). It is also interesting to study the other extreme of 
the hnear (extreme prolate) distribution of (1.2): 

f{p,x) = F{p,)6^'\p^). (CI) 

Our standard assumption of parity symmetry is here that F{—pz) = F{pz). We will spe- 
cialize to the ultra-relativistic case, where we will find that the distribution (CI) is always 
associated with electric instabihties (in the sense of Condition 2 of Sec. II B) . 

Consider directions of q that do not he in the xy plane, i.e. 9 ^ 7r/2. Prom (2.18), we 
get the self-energy 

J.^3^ X 2 f fiP) 

U'{u;,q) = e / — — 
J p P 

Naively, one may use the S function in the distribution (CI) to do all but the p^ integral. 
Summing over the cases Pz < and Pz > (corresponding to = ±1 for this distribution), 

W^{u;,q)^^,mlYl 

± 

with moo naively given by (2.14) as 




The preceding is slightly naive if F(0) 7^ 0, since then the integral (C4) has a logarith- 
mic, small Pz divergence. In physical situations, the 5^'^\p±) in (CI) would have a small 



ctp> -\- q>p^ 



+ 



(-77^ + l)p'p> 



r]+p-q-ie (-rj + p-q-ieY 



(C2) 



+ 



[-7]^ + 1)(5'^5-'^ 



77±cos^ — ie {—ri±cos9 — iey 



(C3) 
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width Apx, and the small pz divergence would cut off when pz ~ Ap^. Then will be 
finite, growing proportional to F(0) ln[(A^)~^] as Ap_L 0, where ~ P±/P and p is a 
characteristic momentum scale of F{pz). If we are only interested in results to leading order 
in this logarithm (which come from ^ Ap), then the equation (C3) is fine.^^ It should 
then be understood that this and future equations in this section are only accurate up to 
relative corrections of order 1/ ln[(A^)~^]. 

As usual, let q lie in the xz plane. The y direction is then an eigen-direction of II'-' , and 
(C3) gives 

Uyy{u;,q)^ml. (C5) 
There is no instability associated with this polarization. The remaining directions give 

n-(a;,g) = ml, (C6a) 

TTxz/ \ 2 sin 6^ COS 6* 

a;,qr = ml- — , C6b 

rj^ — cos^ 



sin^ 6'(?7^ + cos^ 

°° (772 - cos2 ey 



n-(u;,g) = ml—^^^———^, (C6c) 



or equivalently 



n-(a;, q) = rj^U^\., q) = m^ -n^ ^^^^ , (C6d) 
n^"(a;, q) = r/nO"(o;, q) = mlrf sin 9 cos 9 ^^'^ ^ '^'^ Tnl^^ > (C6e) 



2 [-1 + vi'^f coiP' 9 + ri'^ sin* 9 



In all cases, we have absorbed the ie prescription into 77 so that r] now stands for r/ + ie. 

We shall look for all solutions to the dispersion relations in a moment. It is worth noting 
first that n"'"'(0, q) = ml/ cos^ 9, so that Condition 1 of Sec. II B for a magnetic instability 
is not satisfied for ^ 7^ 7r/2. This makes qualitative sense because, for q not in the xy plane, 
the distribution (CI) does not give any trapped particles. [The singularity of n"''^(0, q) for 
9 = 71/2 is a reflection of the zero width of the transverse 6 function in the distribution (CI) 
and will be cut off by a small width Ap±. From Condition 1-b, we know there must be a 
magnetic instabihty associated with 9 — 7r/2.] 

In contrast, n°°(0,q) = m^tan^^ and n°"'(0, q) = do indicate the existence of an 
electric instability by Condition 2 when 

q < gmax = moo tan 9. (C7) 

Note that there is no instability for sin^^ = 0, which is when q points in the same direction 
as the particle motion. As discussed in Sec. II C, the lack of an electric instability in this 



The non-relativistic case is very different. Doing a similar computation of n°''(0, q) in the non-relativistic 
case, one finds a hnear rather than logarithmic small Pz divergence when F{0) ^0. As a result, the 
(5-function approximation in the distribution (CI) is never valid for calculating n°°(0,q) in the non- 
relativistic case when F{0) ^ 0. 
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case is a reflection of the fact that the speed of ultra-relativistic particles cannot change 
from V = 1. 

The nice thing about the distribution (CI) is that we can in fact completely solve the 
dispersion relation in simple closed form. The yy dispersion relation —u!'^+q^+Wy{u!, q) = 0, 
with (C5), just gives the two propagating solutions 



2 2,2 

OJ =q +m^. 



For polarizations in the xz plane, we turn to the xz subspace of the form (Al) of the 
dispersion relation. There will be a solution when det D (taken over this subspace) vanishes. 
Using the results (C6) for n, this condition D'^'^D^^ - {D'^'^Y = [or equivalently D^^D^^ - 
(^9^)2 = 0] can be simphfied to 

g'(r/' - cos' e^irj^Q^ - - 1) + sin" ^ = 0, (C9) 

where we use the dimcnsionlcss variable Q = q/m^. This condition can be factorizcd into 

(^2 _ ^2 _ g-2) JQ2^4 _ g2^2^^ ^ ^2^ ^ ^^2g2 _ ^2^] ^ q (^10) 

where c = cos^ and s = sin 9. There are three corresponding solutions for cu^, in addition 
to (C8), which are 



UJ 



q cos 9 + m^, 



and 



g^(l + cos^ &) ± qsm.Q\j q^ sin^ 9 + 4m^ 



(Cll) 



(C12) 



The lower sign in the last solution gives an instability (cj^ < 0) when q < gmax, as given by 
(C7). This is the only instability of this distribution iox 9 ^ /2. 

APPENDIX D: DISPERSION RELATION FOR PLANAR DISTRIBUTION 



In this appendix, we briefly describe the derivation of the dispersion relation (3.23) used 
to study electric instabilities for the planar distribution (1.1). We look for solutions to the 
matrix dispersion relation (Al) for polarizations in the xz plane by requiring the determinant 
of D in that subspace to vanish: 

= L>^^L>^^ - {D'^'f = (-a;2 ^ ^2 ^ n^^)(-u;' + + XT') - {-q^q, + ^^^)^ (Dl) 
or equivalently 



= D^^D^^ - (D^^f 



-UJ 



+ + g2 + n^") - (n^")2. 



(D2) 



Starting from (3.6) and performing the integral over the angle of p in the xy plane, one 
obtains 



c 
s 



1 - 



cMl 



-1/2 



- 1 



-3/2 



-1/2 



I 2 2 

+ c - 



(D3) 
(D4) 

(D5) 
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where s = sin^ and c = 0086*. Combining these expressions gives a result for (Dl) that is 
times the right-hand side of (3.23). 
In the hmit sin^ — > 0, the matrix D becomes diagonal with 

^ml- rjV, (D6) 
^ q\l -rj') + '^(^-^^+iy (D7) 

In this limit, the unstable mode, given by (3.24), corresponds to E polarized in the x 
direction. Note, in contrast, that the Penrose condition that predicted the existence of this 
unstable mode via Condition 2 was associated with the q eigen-dircction of W^{ie,q) (see 
Appendix A 2), which for sin 6* — > is the z direction. This just reflects the fact that the 
eigen-directions can change as one varies cu from ie to the location i'j of the actual unstable 
solution. 
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Erratum: QCD plasma instabilities and bottom-up thermalization 

[JHEP 08 (2003) 002] 

Peter Arnold, Jonathan Lenaghan, and Guy D. Moore 

In the last sentence of the paragraph containing (3.11), g ~ "^^(i ^i^^ ^)^'^^ should instead 
read q ~ moo[2/(3 sin^ 6)]^/^. 

Before eq. (4.1), there is a misleadingly general statement about when one can ignore 
non-linear interactions induced by hard particles. They can be ignored for gauge fields 
which only depend on a single coordinate, such ets A — A{z) as considered subsequently, but 
not in general. See, for example, Ref. [36], as well as Ref. [40] below. 

Eq. (4.1) should have an overall minus sign on the right-hand side. In eq. (4.2), the e''"'^ 
in the term should be e"^'^. Furthermore, the left-hand side should of eq. (4.2) should have 
an additional term +g€^^'^A^{d'^A'^^ — d^A'^^). The last error is not merely typographic but 
was due to simultaneous cosmic ray events in the brains of the three authors. As a result, 
the solution presented in section 4.1 is not in fact a solution to the equations of motion. 
This does not affect the rest of the paper. Actual static (but unstable) solutions analogous 
to the non-linear solutions of Ref. [36] may be found in Ref. [40] . 

The last paragraph in section 4.1 is misleadingly general. It applies only to cases where 
the typical momenta of unstable modes are of order moo- For the approximately planar 
hard-particle distributions considered in most of the paper, the typical unstable momentum 
scale is instead of order gmax ^ ^oo, where gmax is given by eq. (3.13). 

[40] P. Arnold and J. Lenaghan, "The abelianization of QCD plasma instabihties," hep- 
ph/0408052. 
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